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Abstract: Spatial environmental processes often exhibit dependence in their large values. In order to 
model such processes their dependence properties must be characterized and quantified. In this paper we 
introduce a measure that evaluates the dependence among extreme observations located in two separated 
regions of locations of M 2 . We compute the range of this new dependence measure, which extends 
the existing A-madogram concept, and compare it with extremal coefficients, finding generalizations of 
the known relations in pairwise approach. Estimators for this measure are introduced and asymptotic 
normality and strong consistency are shown. An application to the annual maxima precipitation in 
Portuguese regions is presented. 

1 Introduction 

Natural models for spatial extremes, as observed in environmental, atmospheric and geological sci- 
ences, are max-stable processes. These processes arise from an infinite-dimensional generalization of 
extreme value theory and date back to de Haan (1984) [JJ, Vatan (1985) |17| and de Haan and Pickands 
(1986) [5], who obtained, among other results, a spectral representation of such processes. 

Max-stable processes can be, for example, good approximations for annual maxima of daily spatial 
rainfall (Smith [J3], Coles |3], Schlather |12j . among others) and therefore have been widely applied to 
real data. 

Briefly, a max-stable process Z = {Z x } xeK d is the limit process of maxima of i.i.d. random fields Y x , 
x E TR d , i = 1, . . . , n. Namely, for suitable a n (x) > and b n (x) E H, 

Zx = lim vr = i^ w - w ? x G Rdj 

n->oo CL n {x) 

provided the limit exists. 

We shall consider d = 2, that is Z = {Z x } xl - R 2 ■ The distribution of (Z X1 , . . . , Z Xk ) is a multivariate 
extreme value (MEV) distribution G and we can assume, without loss of generality, that the margins of 
Z have a unit Frechet distribution, F(x) = exp(— x _1 ), x > (Resnick [10|). The distribution G can 
then be defined by 

G(zi,. ..,z k )= exp(-V z (zi, . . -,z k )), Zi E R + , 
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where V z denotes the dependence function of the MEV distribution G, which is homogeneous of order 
-1, i.e., V x (azi, ... , az k ) = a~ 1 V r z (zi, ... ,z k ), z { G R + , a > 0. 

The dependence function captures the multivariate dependence structure and the scalar V z (l, . . . , 1) 
defines the extremal coefficient considered in Schlather an Tawn |13] which measures the extremal depen- 
dence between the variables indexed in the region z. This coefficient varies between 1 and k depending on 
the degree of dependence among the k variables. These measures of dependence have gained great impor- 
tance since quantifying dependence between extreme events occurring at several locations of a random 
field is a fundamental issue in applied spatial extreme value analysis. 

Cooley et al. (2006) showed that the bivariate extremal coefficient can be directly estimated from the 
madogram (that represents a first order variogram), obtaining in this way a connection between extreme 
value theory and the field of geostatistics. A measure of the full pairwise extremal dependence function 
is given by the A-madogram defined in Naveau et al. [9j as 



where F denotes the marginal distribution of Z. 

Although the A-madogram fully characterizes the pairwise extremal dependence it does not enable 
the analysis of dependence between maxima over two disjoint regions of locations. The importance of 
characterizing dependence between extremes occurring at two disjoint regions has been recognized by 
hydrologists who have grouped data into regions based on geographical or catchment characteristics. It 
is clear that, for example, the weather in mountain regions usually affects the weather in the surrounding 
regions and therefore a measure able to capture such regional dependence is essential. 

In this paper we propose a measure that enables the analysis of dependence between maxima over 
two disjoint regions of locations x — -fjci, . . . , x k } and y = {y 1} . . . , y s } and therefore generalizes the 
A-madogram. This measure, here called generalized madogram, is introduced in Section 2 and some of 
its main properties are presented, namely its relation with the dependence function of the MEV G. In 
Section 3 we present estimators for the generalized madogram and derive the respective properties of 
strong consistency and asymptotic normality. The performance of the proposed estimators is analyzed 
in Section 4 with a max-stable M4 random field. Finally, Section 5 illustrates our approach through an 
application to precipitation data from Portugal. 

2 Generalized madogram and dependence of spatial extreme events 

The measure that enables the analysis of dependence between maxima over two disjoint regions of 
locations x = {x±, . . . , x k } and y = {y 1: . . . , y s } and generalizes the A-madogram introduced in Naveau 
et al. [9] is defined as follows. 

Definition 2.1 Let Z = {Z :r } 3 , gR 2 be a max-stable random field with unit Frechet margins and x = 
{x±, . . . , x k } and y = . . . , y s } two disjoint regions of ]R 2 . The generalized madogram is defined as 




(1) 




a > 0, > 



(2) 



where M(x) = V tl Z m and M (^) = Yj=i Z ; 
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Remark 1 The following equalities hold for the generalized madogram 

k s 



z/^(x,y) = & 



1 



-E 



V F%Z Xi )-\l F?(Z yj 

3=1 



i=l 



1=1 



a 



This representation of the generalized madogram, i^ a '^(x,y), will motivate the definition of natural esti- 
mators, in Section 3, for this coefficient. 

Remark 2 When we take /3 = 1 — a, a 6 (0, 1), and k = s = 1 in (j2J), we obtain ([I]). 

The following proposition states that z/ Q, ^(x, y) provides dependence information between the regions 
x and y through the dependence function of the MEV distribution G. This result generalizes Proposition 
1. in Naveaux et al. [9\. 

Proposition 2.1 For any max-stable random field with unit Frechet margins and for each pair of disjoint 
regions of locations x = {xi, . . . , x k } and y = {y 1 , . . . , y s } in K 2 , we have 

V Xt y(a,.. .,«,/?,... ,/3) 



^(x,y) 



1 + V Xi y(a, . . . ,a,(3, . . . ,/3) 



c(a,P) 



with 



where 



and 



2 V« + F x (i,...,i) /3 + y y (i,...,i); ' 

G x ,y(zx,...,z k+S ) = P ( If] {Z Vi < Zi \ X f] \f] {Z Vi < z k+i }X) , zt € M+. 



J=i 



.i=i 



Proof: To obtain the result, we start by transforming the definition of i^ a, ^(x, y) through the relation 
\a — b\ = 2(o V b) — (a + 6), and then take into account that 

E (F* (M(x)) V (M(y))) = E (f V ^M)) . 



It holds 



ft 



/3 



P(M(x) < «n,M(y) < /3u) 

G Xi y(o;it, . . . , cm, f3u, . . . , /3m) 

exp {— u~ 1 V x ,y(a, a, /?,..., /3)} , -u > 0. 
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Hence, 



+0 ° d 

F(u) exp (— V K>y (au, . . . , au, f3u, . . . , (3u)) — (— V X{ y(a.u, . . . , era, f3u, . . . , f3u)) 

r+oo 

= / exp - u _1 y Xiy (a, ...,a,P,...,{})) u~ 2 V Xi y(a, ...,a,P,..., P)du 

Jo 

V x ,y{a, ...,a,P,...,P) 
l + V x , y (a,...,a,P,...,P) 

Using similar arguments, we obtain E(F a (M(x))) = ^fe 112 ^ t = ^vTTij an< ^ 



Remark 3 For each a, {3 > the coefficient c{a, /3) considers the dependence intra each of the regions 
x and y through the extremal coefficients of vectors with margins Z Xl , . . . , Z Xk and Z yi , . . . , Z Vs . When 
we consider c(a,P) constant, the dependence between x and y is stronger for lower values of z/*'^(x, y), 
corresponding to lower values of V x>y (a, . . . ,a, P, . . . , P). 

In the following proposition we establish some properties of the generalized madogram. 

Proposition 2.2 Let x = {x±, . . . , Xk} and y = . . . , y s } be disjoint regions of H 2 . We have, for 
each a, [3 £ IR + , 

1. 0< v a >P(x,y) < |; 

2 - ^' a (x,y) = - \ (sfe + 5^) , where e xUy = U x , y (l, . . . , 1). 



Proof: The first statement results from the definition of the generalized madogram and the second 
follows from the definition of the extremal coefficient e. □ 



Remark 4 The function z^ a ' a (x, y) can also be related with the dependence coefficients considered in 
Ferreira [7j as follows: 

a,at \ e y £ i( x .y) / v (e y + e x )e 2 (x, y) 

v ' (x, y) = — f - c(a, a) = — f- ■ r— - c(a, a), a > 0, 

a + e y ei(x,y) a + (e y + e x )e 2 (x, y) 

where ei(x, y) = xUy and 62 (x, y) = — xUy . These coefficients evaluate the strength of dependence 
e y e x + e y 

between the events {M(x) < u} and {M(y) < u} . 
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Remark 5 If the variables M(x) and M(y) are independent then 

e x + e y 1 / e x e y 



^ Q (x,y) = * y - £ — 2- + 



a + e x + e y 2\a + e x a + e y 

whereas if they are totally dependent 

a a / \ V 6y 1 / £ x £y 

z/* >Q! (x,y) = — -f — — + 

3 Estimating the generalized madogram 

Proposition 2.1 relates the generalized madogram with well known dependence measures. Immediate 
estimators for the generalized madogram can then be obtained through the estimators of those measures, 
which have already been studied in the literature. For a survey we refer to Krajina jS], Beirlant pQ and 
Schlather and Tawn [13] . 

In this section we present a natural non-parametric estimator for the generalized madogram based on 
sample means. 

Let (Zxi , • • • , Ex}) an d (Zy± , ■ ■ ■ , Zy})i t = 1) ■ ■ ■ )T, be independent replications of (Z Xl , . . . , Z Xk ) 

and (Z yi ,...,Z y J, respectively. Hence {M t (x) = \J k i=1 Z$, t = 1,...,T} and {M t (y) = Vi=i z % \ 

t = 1, . . . , T} are random samples of M(x) and M(y), respectively. 

In the case of known maginal distribution F, which becomes unit Frechet by transformation —-, — } lr7 ; 

o > J log F(Z Xi ) 

for a^EK 2 , the estimator for the generalized madogram is given by 

1 T 

^( x .y) = ^E™( x ))- F/? ( M »W)l' «>o, /3>o. 

1=1 

This estimator is unbiased and converges in distribution to a Gaussian distribution, as stated in the 
following proposition. 

Proposition 3.1 (Asymptotic normality and strong consistency under known marginal distribution F) 
We have 

VT(^(x,y)-^(x,y)) ^ 

(7 

w/jere a 2 = i 7 »^( x ,y) - (z/^(x,y)) 2 and 7^(x,y) = \E [(F a (M(x)) - F^(M(y))) 2 ] . Moreover, 
z/^^x, y) converges almost surely to z/*'^(x, y). 

Proof: Let Yi, . . . , Y T be independent copies of Y = \ |F a (M(x)) - F /3 (M(y))|. We have that 

ay 

where M y = \E |F a (M(x)) - i^(M(y)) | = i/°^(x,y) and 4 = ^W) " (^(x, y)) 2 . 

The strong consistency of ^"'^(x, y) follows since the sample mean converges almost surely to the 
mean value. □ 
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When the distribution of each Z Xi , F Xi , is unknown we take the empirical Frechet normalization of 

the variables, i.e., uj* = ^ — ttt— - 1 • ' ''" ' ; 

find the following modification of the above estimator 



^ — TjT— , where F Xi (Zx/) is the empirical distribution function. We then 

^g(F Xi (Z^.')) 



T 

v (x,y) = ^E 



3=1 



\i=l 



where 



1 T 

-Y 



2T 



t=i 

T 



(*)' 



a 



2T ^ 

t=i 



V K (4>) - V *S «) 



r^-{zw<u}- 
t=i 



a > 0, /3 > 0, 



This estimator, z? ' (x, y) is strongly consistent since it holds 



T 



t=l 



i=l 1 = 1 



(*)■ 



2 



Vf-(<)-v^) 



i=i 



< 



T 

-E 



t=l 



i=l 7 = 1 



T 

--E 



i=l 



i=i i=i 



T 

9 T ^ 



t=l 



yF»{Z®)-\jFf>{Z®) 
i=i j=i 

The second term converges almost surely to zero by the strong Law of Large Numbers. 
In what concerns the first term we have 



2 



V F°{Z$) - V FP{Z$) 

i=l j=l 



T 

--Y 

9 T 



t=l 



Mn(z$)-\iF% j{ zi 

i=l 7=1 



T 

2T^ 



t=i 



V^ Q (^)-V^(4?) 



i=i 



< -E 



t=l 



i=l 7 = 1 



Vf a (^)-V^(4! 



i=i 



< 



T 

9 T ^ 



t=l 



V *5 0$) - V - V ^(^) + V ^ 

1=1 1 = 1 j=l 1 = 1 



(t)< 



< 



< 



1 1 

2T 



1 1 
2T 



E 

t=i 

T 

E 

t=i 



i=l j=l 



E 0$) - ^(^)| + E |^(^) - ^ 
i=i j=i 
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which converges almost surely to zero. This follows from the strong consistency of the empirical distri- 
bution function and the fact that y Y2i=i l-^il converges almost surely to where Xj, i = 1, . . . ,T, 
are i.i.d. random variables with -ElXil < oo. 



The asymptotic normality of the estimator is obtained by taking J(x\, . . . , x k+s ) = | Vj=i x ? ~ Vi=fc+i x 
in the following theorem stated in Fermanian et al. [6]. 



Theorem 1 Let (Z Xl , . . . , Z Xk , Z y , . . . , Z y ) be a random vector with d.f. H and continuous marginal 
d.f. 's F X1 , . . . , F Xk , F y , . . . , F y and let the copula Ch have continuous partial derivatives. Assume that 
J is of bounded variation, continuous from above and with discontinuities of the first kind. Then 

1 T 

—f= { J ( F xi( z $)i ■ ■ ■ i F x k (Zxk)> F yi( z Vi)i ■ ■ ■ > F y s ( z ys)) 

- E(J (F Xl (Z X % l ),..., F Xh (Z X % l ) , Fy x (Zy 1 } ),-.., Fy g (Zy^ ) | 

G(ui, . . . ,u k+s )dJ(u 1 , . . .,u k+s ), 

[o,i] fc+s 

where the limiting process and G are Gaussian with zero mean. 

In the following section we will conduct a simulation study of an M4 random field to assess the 
performance of the estimator ^ Xj y 

4 An M4 random field 

It is well known that the class of max-stable processes called multivariate maxima of moving maxima 
processes or simply M4 processes, introduced by Smith and Weissman |15| . is particularly well adapted 
to modeling the extreme bahaviour of several time series. 

To illustrate the computation of the generalized madogram given in ([2]) we will now define an M4 
random field. 

Lets consider that the distribution of (Z Xl , . . . , Z Xp ) is characterized by the copula 



hoc +oo 



C(u Xl ,...,u Xp ) = H H /\ u a J™, u Xi G [0, 1], i = 1, . . . ,p, (3) 

1=1 m=-oo a:g{a;i,...,Xp} 



-oo +oo 



where, for each x G Z 2 , {ai ma; }i>i,mez are non-negative constants such that a,i mx = 1. This 

Z=l m=— oo 

random field Z is max-stable, since, for each t > 0, the copula (J3]) satisfies 



C f (u xi ,. . . ) = C(u 



.(',, ' v S "■(•!> • • • ) U Xp)l 



for any locations xi, . . . , x v . 

As the M4 process considered in Smith and Weissman [15], we can consider that for each location x, 
Z x is a moving maxima of variables Xi tTl , i.e., 

Z x = max max ai mx X L1 _ m , x G Z 2 , (4) 

Z>1 — oo<m<+oo 
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where {-2Q,n}z>i,nez is a family of independent unit Frechet random variables. The dependence structure 
of (Z X1 , . . . , Z Xp ) is regulated by the signatures patterns ai mx and is given by ([3]). 
For each pair of regions x = {x\, . . . , x k } and y = {y k +i, ■ ■ ■ , Vk+s} we have 

V Xj y(zi,...,z k ,z k+ i,...,z k+s ) = -luC(e'^ 1 ,...,e~ z k+s) 

+00 +00 k+s 

= ^2 V ^"-imx,, z% G E,, i = 1, . . . ,k + s, 

1=1 m=—oo i=l 

and consequently, for a > and /3 > we obtain 

+00 +00 / k k+s \ 

Vx, y (a, ... ,a,(3, ... ,f3) = ^ ^ I \/ a -1 a/ mxi V \f $-~ x a lmx . J . 

i=l m=-oo \i=l j=fc+l / 

To illustrate the computation of the generalized madogram we shall consider, in what follows, examples 
with a finite number of signature patterns (1 < I < L) and a finite range of sequencial dependencies 
(Mi < m < M 2 ). 

Example 4.1 Lets consider that for each location x 6 1? with even coordinates we have 
o-iix = cbvzx = \ and otherwise auj. = 4 = 1 — ai2a:- The values of (an x , ai2 X ) determine the mov- 
ing pattern or signature pattern of the random field, which in this case corresponds to one pattern 
(L = l). 




Figure 1: Simulation of the M4 random field with L = 1, 1 < m < 2 and for constants defined in Example 4.1. 



For the disjoint regions of locations x = {(2, 1), (2, 2)} and y = {(3, 3), (3, 4)} we have 

V x , y (a, a, ft p) = i(2a- 1 V /T 1 ) + ^(a' 1 V /T 1 ) 



and therefore, the generalized madogram in this pair of locations is given by 



Example 4.2 Lets now assume that for each location x = £ Z 2 , an x =4 = 1 — ay^x if i < j and 

ailx = | = 1- a-i2x if i > j- 

As in the previous example the M4 random field generated by these sequences has a single signature 
pattern. Considering now two disjoint regions of locations with different size, x = {(1, 1)} and y = 
{(3, 2), (3, 3), (4, 3)}, we obtain 

V x>y (a, (3, p, P) = \{a~ l V 3/T 1 ) + -^a' 1 V /T 1 ) 
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and consequently 



(x,y) 



Ua- 1 V3/T 1 ) + vr 1 ) 



1+ (a- 1 V3H + 



a 



V/?- 1 ) 



+ 



a + 1 /3 + 



a > 0, p > 0. 




j i i 

Figure 2: Simulation of the M4 random field with L = 1, 1 < m < 2 and for constants defined in Example 4.2. 



Example 4.3 As stated in Zhang and Smith |18j . in a real data generating process it is unrealistic to 
assume that a single signature pattern would be sufficient to describe the shape of the process every time 
it exceeds some high threshold. Hence, we shall now consider one example with two signature patterns 
(L = 2). 

Lets assume that for each location x = we have a\\ x = a\ 2x = «i3a; = t^, 

a2ix = a-22x = a23x = i if both coordinates are odd and au x = jg, a\ 2x = \, a\3x = g, 
a-2ix = a 2 2x = a 2 3x = § otherwise. Now the values of (a nx , a Ux , ai 3x ) and (a 21x , a 2 2x, a 23x ) define 
the two signature patterns of the random field. 




Figure 3: Simulation of the M4 random field with L = 2, 1 < to < 3 and for constants defined in Example 4.3. 



For two disjoint regions x = {(2, 1), (2, 2)} and y = {(2, 3), (3, 3)} we now have 

V x , y (a, a, (3, (3) = (c^ 1 ! V /T 1 ^) + ± (a' 1 V /T 1 ) + \ (a" 1 V /T 1 ) + (a _1 | V /T 1 ^) 



and consequently 

z/^(x,y) 



18 v 12 ,/ t 9 "i" 6 " l "^3 v 4 
x ^^18 v 12y^ 9 "•" 6 ' ^ 3 v 4 

if l 10 



These examples will be used in the following simulation studies to assess the performance of the esti- 
mator given in Q. Figures 4, 5 and 6 show the simulation results obtained by generating 50 replications 
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of 100 independently and identically distributed max-stable M4 random fields in the three situations 
previously presented, with a and (3 taking values in {k x 0.2 : k = 1, . . . , 100}. 

As we can see from the values of the mean square error (Figures 4, 5 and 6) the estimates obtained 
from our estimator v (x, y) are quite close to the true values of the generalized madogram. 




20 



Figure 4: The true values of the generalized madogram (z/*'^(x, y)), the estimated values (D (x, y)) and the 
estimated mean squared error (MSE) for Example 4.1 (x = {(2, 1), (2,2)}, y = {(3,3), (3,4)}). 



^"'^(x.y) d a ' (ic.,y) MSE 




Figure 5: The true values of the generalized madogram (z/*>^(x, y)), the estimated values (P (x, y)) and the 
estimated mean squared error (MSE) for Example 4.2 (x = {(1, 1)}, y = {(3, 2), (3, 3), (4, 3)}). 



j/"^(x,y) 5 Q,/3 (x, y) MSE 




Figure 6: The true values of the generalized madogram (z/*'^(x, y)), the estimated values (D ' (x,y)) and the 
estimated mean squared error (MSE) for Example 4.3 (x = {(2, 1), (2, 2)}, y = {(2, 3), (3, 3)}). 



5 Application to precipitation data 

We now consider an application of the proposed estimator of the generalized madogram to annual 
maxima values of daily maxima rainfall in different topographic regions of Portugal. 

In Portugal there are topographic differences between north and south, whereas in the north there 
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are several mountains, in the south there are mainly plateaus and plains. The Central Cordillera, formed 
by the mountains of Sintra, Montejunto and Estrela, divides south and north and influences climate in 
Portugal, namely precipitation. The majority of the precipitation in Portugal comes from North- West 
and it is more abundant in the north than in the south due to the Central Cordillera that creates a 
physical barrier for precipitation. 

To study the dependence between extreme precipitation occurring in this mountain area and in 
surrounding areas, we considered precipitation data that consist of annual maxima of daily maxima 
precipitation recorded over 38 years (1944-1981), in 5 Portuguese stations, obtained from the Portuguese 
National System of Water Resources (http://snirh.pt). In Figure [7] we can view the location of the 5 
stations considered. We remark that the stations "Lagoa Comprida" and "Fajao" are located in North- 
West part of "Serra da Estrela" , the highest mountain in continental Portugal and part of the Central 
Cordillera. 



Ponte 
Lime 

Varzitn Q Braga 
Porto O * p sr?<l«s 




Pomtial Coimta V 

o Castelo 
Portugal Bra n« 

o 

Caret! 



CascajsQ O Lisijoa 




(223368.165.390619.802) 
Caktas da Felgueira 



(282807.355,35531 7.B08) 



(279821.567,330915.754) 
Barragem Cabccc Monte o 



200000 220000 240000 260000 280000 300000 
x(m) 



Figure 7: The locations of the stations where precipitation data were collected, obtained from Portuguese 
National System of Water Resources (left) and their representation in Lambert coordinates (right). 



Since the data are maxima over a long period of time, we assumed that they are independent over 
the years in each location. We also assumed that the random field is max-stable with unknown marginal 
distributions so data were previously transformed at each site so that they have a standard Frechet 
distribution. 

In Figure 8 we picture the estimated values of the generalized madogram, for several values of a and /3, 
when considering the mountain region x = {Fajao, L. Comprida} and different topographic surrounding 
regions, that are either to north or to the south of this region. 

The estimated values for a = X/k and j3 = (1 — A)/s, A £ (0, 1), where k and s are the number of 
locations in region x and y, respectively, are presented in Figure 9, for the three situations considered in 
Figure 8. 

Since lower values for z/*'^(x, y) indicate strong dependence, the results presented in Figure 9 suggest 
a stronger dependence between the mountain region x = {Fajao, L. Comprida} and the north region 
y = {C. Felgueira}. This is in accordance with the previously stated that the Central Cordillera creates 
a physical barrier for precipitation in Portugal. 

It would be interesting to further investigate this dependence with other regions but the lack of 
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x — {L.Comprida,Fajao} y — {B.C.M.} x — {L.Comprida,Fajao} y — {Pcnamacor} x — {L.Comprida,Fajao} y — {C.Fclgucira} 




0.0 0.2 0.4 0.6 0.8 1.0 

X 

Figure 9: Generalized madogram id ' (x, y)) estimates obtained for each pair of disjoint regions when a = -| and 
P = 1- A with A e (0,1). 

available data restricts the possible regions to study. 

All the simulations presented in this paper were done in R statistical computing program 
( |http:/ /cran.rproject) org/). We remark that several packages on Extreme Value analysis have been 
recently introduced into R, but more recently Ribatet [TT] added to R the package SpatialExtremes that 
provides functions to analyze and fit max-stable processes to spatial extremes. 
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